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ABSTRACT 

We exhibit sample parameter sets for the SO(10) NMSGUT[U [2] superpotential, 
and associated Susy soft spectra that match, at low energies, with SM fermion(includin^ 
neutrino) precision data and are consistent with Unification constraints including 
Mb-l ~ Mx, besides yielding Right handed neutrino masses in the Leptogenisis rel- 
evant range 10 8 — 10 13 GeV. Matching the precison data requires lowering strange and 
down quark masses above Ms via large tan/? driven supersymmetric threshold cor- 
rections : which requires characteristic soft Susy breaking spectra. Thus our solutions 
evade the numerical 'no-go' objection raised by Grimus and Lavoura and exactly real- 
ize our proposal[3] that tiny 126 couplings allow realistic neutrino masses even though 
M B _ L ~ 10 16 GeV : as is generic in single scale SO(10) Susy GUTs. The /i and trilin- 
ear soft A t parameters are positive and negative respectively and one or more third 
generation sfermions tend to be heavier than the first two generations: this marks out 
a distinctive NMSGUT experimental signature in the sfermion spectrum potentially 
discoverable at LHC. We estimate the d — 5 RRRR decay rate for the dominant 
p — > K + u T channel. It indicates that d — 5 baryon decay rates will further winnow 
candidate parameter sets and strongly limit sfermion mixing angles. The NMSGUT 
is thus reaching the stage where it's parameters are subject to enough constraints to 
overdetermine them and thus subject it to stringent falsification tests: given sfermion 
data from the LHC and from on going flavour and B violation experiments. 



1 Introduction 

Grand Unification has long retained its attractive status as the most audacious ex- 
trapolation of the structural logic of the Standard Model which is not so hyperbolic 
as to have lost every possibility of ever satisfying common and well established de- 
mands of the scientific method (such as falsiflability or even verifiability) . It enjoys 
the support of all of the few clear clues about beyond Standard Model(SM) physics 
that have so far emerged. Yet only infrequently, if at all, have there been attempts 
to fully specify a Grand Unified model (GUM) at a level of detail approaching the 
SM, so as to be able to confront the fully specified model with all the experimental 
data in hand and to extract the consequences thereof without the aid of egregious 
assumptions. The commonly accepted necessity of invoking Supersymmetry in order 
to ensure the stability of the parametrization of physics in the extreme Ultravio- 
let has drawn an additional veil, woven of the plethora of unknown supersymmetry 
breaking parameters, over the portal to the Grand Desert at the scale(Ms) where 
supersymmetry is (softly) broken. As we enter the epoch of the 'discovery machine' 
LHC we still lack well defined candidate patterns of Ms ~ 1 TeV scale physics which 
are well correlated with any specific supersymmetric GUM. In this letter we begin 
to delineate emergent patterns associated with the so called minimal renormalizable 
Supersymmetric SO(10) models [11 15| IB"! \7\ l9l l8| ITTj [3] which are arguably the best mo- 
tivated and well defined of supersymmetric GUMs. By numerically investigating to 
what extent the mass formulae and unification constraints, written directly in terms 
of the fundamental GUM parameters'^, can accommodate the known fermion mass 
data, the constraints of unification and those from limits on exotic processes we arrive 
at tentative information about soft supersymmetry breaking parameters required to 
prevent the NMSGUT from failing to fit the available fermion mass data, unification 
constraints and limits on exotic processes. 

The question whether the renormalizable Supersymmetric SO(10) GUT based on 
the Higgs set 210 © 10 © 120 © 126 © 126 (i.e the so called Next or New Minimal 
Supersymmetric GUT NMSGUT pQ ( although a more accurate name might be the 
Full Minimal Supersymmetric GUM(FMSGUM) !) is capable of encompassing fermion 
mass-mixing data has received much attention lately [TH [31 [121 131 El EH El E]. In 
particular the demonstration that generically feasible neutrino mass patterns could 
not actually be realized[TTl [3] in the classic MSGUT[U [5] ( in the context of the 
Babu-Mohapatra[T7j program of restricting the FM Higgs irreps to just the 10, 126) 
dramatically highlighted the possibility and necessity of using GUM specific mass 
formulae formulae. In particular the generic sub-dominance [TT| [3] of the Type II 
seesaw in the MSGUT suggested^ [121 E] that by using very small 126 couplings to 
lower the right handed neutrino masses one might enhance the Type I Seesaw masses 
which were otherwise generically ~ v^/Mx < 1 meV and thus too small to account 
for the measured neutrino oscillation data. This would obviously be possible only if 



the CKM structure were generated by the third possible Higgs representation i.e the 
120. However this idea almost immediately faced an apparently strong challenge |18j 
which showed (numerically) that the fit to charged fermion data using only 10 and 
120 representations always failed badly to account for the down quark masses which 
always come out much too small if the 3rd generation was fitted correctly. It is worth 
mentioning here that till recently these discussions -in the 'MSGUT community' 
that is- proceeded under the somewhat sanguine assumption that the appropriate 
uncertainties for the fitting exercise could be estimated by ignoring all new threshold 
effects right up to the GUT scale and hence were essentially the extrapolated errors 
at the low scale : thus for example the lepton masses were sought to be fitted to 
one part in a million or better [38 j. In fact however it has been known [T9| 120] from 
the early days of large tan/? b — r — t unification that the T 3L = —1/2 charged 
fermions suffer large threshold corrections and the relevance of these corrections for 
estimating the uncertainties in the GUT scale fermion masses has recently been re- 
emphasized [221 123] . 

However this was followed by other work in which we found [T2| [T3] that our 
scenario implied somewhat recognizable simulacra of the fermion data and the same 
authors [TU [15] even showed that generic fits with small 126 couplings but much larger 
10, 120— plet ones could in fact accurately fit all the fermion data available to date 
: thus making their own earlier and quite cogent generic arguments somewhat moot. 
However the question of whether the specific structure of the NMSGUT might still, 
as in the MSGUT case, defeat the demonstrated generic possibility remained open. 
In[2] we showed that in order to obtain fits that were compatible with the NMSGUT 
it would be necessary to search for fits using the detailed NMSGUT specific formulae 
developed by us in [1] since the the generic parameters used in [T^ [T5] were in fact 
subject to highly non-trivial constraints (due to their origin in the NMSGUT) that 
were in practice unverifiable since they could not be inverted. Two recent papers [22] 
123] re-emphasized that down type fermions could suffer very large corrections to their 
masses at the Supersymmetry breaking scale Ms and that the formulae for computing 
these effects to all orders in tan /3 were long known to experts if not so well appreciated 
by the MSGUT community [201 HI EI]- This reminder [23J came at a particularly 
opportune time for our numerical investigation of the fitting problem in the NMSGUT 
which, while always tumbling naturally towards ultra low 126 coupling values, was 
constrained, once it had satisfied the strong demand for accuracy in the top quark 
and charged lepton masses, to settle for very low d and s quark yukawa couplings at 
the GUT scale. This was indeed the very behaviour noticed-and judged pathological- 
in [IE]. However the new emphasis on large tan/? uncertainties in these very masses 
naturally suggested to us that the NMSGUT (which like other £0(10) GUTs prefers 
large tan/3 to naturally explain third generation charged fermion yukawas as nearly 
equal due to their domination by the 33 component of the 10— plet alone, as is 
possible -and highly persuasive- precisely because all of them lie in a single 16- 



plet) could readily avail of this crutch that seemed almost designed to aid its frailty. 
Taking this new found liberative message at face value, and the difficulties of the 
NMSGUT with d, s quark masses as a signal from the model of its character and 
inner necessity, we have used the diagonal threshold correction formulae to match the 
down and charged fermion masses to the rather low values found by our "downhill 
simplex" based search and fitting program (after running them down from Mx using 
the full two loop equations of [32] to Ms); we also include running of the Weinberg 
operator [33] [31]. 

The analytic solution of the fitting equations, is practically inconceivable in view 
of their staggering nonlinearity[U [2] even in the case of generic fits, although a 
perturbative approach [121 02] may now well work(see Section 3). Thus even the 
generic demonstration of feasibility [HI [15] had recourse to a nonlinear numerical 
algorithm(the 'downhill simplex' method). As explained in [H[2] the fitting problem 
becomes much more complex when formulated in terms of the actual parameters of 
the underlying GUM theory. Thus all treatments of this problem must resort to nu- 
merical analysis at some stage. Indeed, saving some un-foreseeable miracle, realistic 
consideration of the fitting problem shows that the entire program of understanding 
fermion masses in a GUM context must finally bow before the intrinsic non-analyticity 
implied by the highly non linear relation between the fundamental parameters and 
those of the low energy theory. We have used numerical algorithms (based on the 
so called downhill simplex method of Nelder and Mead [25]) which incorporate our 
explicit (but very complicated) fermion mass formulae [1] in terms of NMSGUT pa- 
rameters and searched extensively for plausible parameter sets for the NMSGUT while 
demanding that the parameters be compatible with both the known fermion data and 
the requirements of unification, besides plausible criteria such as perturbativity. 

In a sense this effort is still premature. Even with the assumption that the param- 
eters of the NMSGUT superpotential are real (save the parameter that is fine tuned 
to keep one pair of doublets light), and thus only 23 in number, there are still 5 more 
parameters than available data {m«, c ,t,md, s ,(),me,/i,T, Amf 2 , Am| 3 , #12,23,13; S CKM , 
#122313} {6i3 MNS < 10° is only an upper limit [26]). The presently available fermion 
data thus comprise only 18 parameters : three leptonic CP phases and the absolute 
value of one neutrino mass remain unknown for the moment. The NMSGUT, on the 
other hand has 23 free parameters. These are however subject to various additional 
structural constraints such as the requirements of unification on as{Mz),Mx and 
of neutrino masses on M^_^(via the Seesaw mechanisms naturally operative in the 
SO(10) GUTs) which are in fact closely related through their dependence on the spon- 
taneous symmetry breaking pattern in the NMSGUT. Additional constraints may be 
expected to emerge when one examines the compatibility of the fundamental param- 
eters with the emergent data. Indeed in the foreseeable future we may have access to 
another Neutrino angle and phase as well as a neutrino mass and also obtain infor- 
mation about the pattern of spartner masses at LHC. Thus the day may not be too 



far when the fitting problem is fully determined or even over determined. At that 
point at least the NMSGUT will mature as a scientific (i.e falsifiable and predictive) 
hypothesis. Thus although the fitting problem is both under determined and ana- 
lytically insoluble the fundamental significance of the problem justifies a numerical 
attack to see if any trends that may inform experimental searches for supersymme- 
try and searches for exotic processes can be extracted when the GUM begins to be 
specified at this unprecedented level of detail. 

From the searches we have so far carried out on PC based computer systems, 
certain persistent and very specific themes have emerged. In particular we find that 
the achievement of adequate neutrino masses and large lepton mixing require that 
the masses of the strange and down quark above the Susy breaking scale M$ must be 
reduced by a factor of 3-4 while raising the b, r masses by 5-30%. As mentioned above, 
such an ad hoc distortion might have seemed unacceptable to some workers in the 
field until the fitting effort in GUTs, and specially in the SO(10) MSGUT [38J, was 
reminded of the large Susy threshold corrections operative at large tan j3 for the down 
type quarks (22J [23]. A fit of the fermion spectrum achievable at the GUT scale in 
the NMSGUT then seems to require a characteristic scenario for Susy breaking mass 
patterns. Since (accurate) unification is the main motivation for Supersymmetry 
itself, the need at this juncture is for GUM driven insight into the possible mass 
patterns of Susy partners that might be unveiled by LHC We offer the preferred Soft 
parameter patterns of the NMSGUT as our 'bet' in the LHC-Susy 'sweepstakes' ! 

In Section 2 we briefly describe the NMSGUT and the fermion fitting formulae 
specific to it. Complete details can be found in pfl [2]. In Section 3 we describe 
our numerical procedure and then give typical fits found by our search algorithms 
and dissect their structure. The required pattern of threshold corrections provides 
indications about the pattern of sfermion masses and other Susy breaking parameters 
and thus constitutes a first (qualitative) prediction of the NMSGUT about potential 
discoveries at the LHC as well as about the nature of the soft supersymmetry breaking 
terms at the GUT scale(after extrapolation back to that scale). In Section 4 we use 
the explicit sample sets of NMSGUT parameter values and Susy soft breaking data 
to estimate typical predictions of a fully specified NMSGUT for the RRRR diagram 
for p — > K + v nucleon decay which is dominant at large tan/5[28j . In Section 5 
we conclude with a discussion about the outlook for discovery as well as the further 
development of the fitting program. 

2 NMSGUT fermion mass formulae 

The NMSGUT [lj is a renormalizable globally supersymmetric SO (10) GUT whose 
Higgs chiral supermultiplets consist of AM(Adjoint Multiplet) type totally antisym- 
metric tensors : 210($yfe/), 126(1]^™), 126(Syjw m )(i,j = 1...10) which break the 
SO (10) symmetry to the MSSM, together with Fermion mass (FM) Higgs 10 (H$) 



and 120(Oijk)- The 126 plays a dual or AM-FM role since it also enables the gen- 
eration of realistic charged fermion and neutrino masses and mixings (via the Type I 
and/or Type II Seesaw mechanisms); three 16-plets ^a_(A = 1, 2, 3) contain the mat- 
ter including the three conjugate neutrinos (u£). The superpotential (see[6], [Tl [HI El H] 
for comprehensive details ) contains the mass parameters 

m : 210 2 ; M : 126 ■ TM; M H : 10 2 ; m : 120 2 (1) 

and trilinear couplings 



A : 210 3 ; T) 

k 

c 



210 • 126 • 126; 7 © 7 : 10 • 210 • (126 © 126) 
10 120 -210 ; p: 120 120 -210 



120 -210 126 ; C : 120 • 210 ■ 126 (2) 

In addition one has two symmetric matrices h^B, Jab of Yukawa couplings of the 
the 10, 126 Higgs multiplets to the 16a-16b matter bilinears and one antisymmetric 
matrix qab for the coupling of the 120 to 16a-16b- One of the complex symmetric 
matrices (which, as in most previous work, we choose to be Hab) can be made real and 
diagonal by a choice of SO(10) flavour basis. Thus initially complex Yukawas contain 
3 real and 9 complex i.e 21 real parameters. If the Yukawas are taken to be real 
however one will have just 12 real yukawas. One out of the rest of the superpotential 
parameters (i.e m, Mh, M, Mo, A, 77, p, k, 7, 7, C, C)> sa Y Mr, can be fixed by the fine 
tuning condition to keep two doublets light. As explained in0, [7J El [HI El Q] the 
fine tuning fixes the Higgs fractions i.e the composition of the massless electroweak 
doublets in terms of the (6 pairs of electroweak) doublet fields in the GUT. A subtle 
point here is that even if the other parameters are taken real the fine tuned Mh 
(which does not itself enter into the low energy lagrangian) will be complex and this 
obscures the issue of whether the resultant CP violation in the low energy effective 
theory should be considered 'spontaneous' or not. Since doublet triplet splitting is 
a common unresolved problem for most GUMs we shall take a pragmatic approach 
and assume all superpotential parameters except the fine tuned Mh to be real. One 
could refer to this induction of phases in the low energy effective superpotential via 
the super heavy vevs and the fine tuning as "effectively spontaneous CP violation "! 

MSGUTs enjoyjElElIB] a crucial simplification : the GUT scale vevs and therefore 
the mass spectrum are all expressible in terms of a single complex parameter x which 
is a solution of the cubic equation 

8x 3 -15x 2 + 14x-3 = -£(l-x) 2 (3) 

where £ = ^. 

If we ascribe the generation of phases in this theory to the high scale symmetry 
breaking ( the fine tuning mechanism that keeps a pair of MSSM doublets light stands 
in for the unknown dynamics that enforces doublet triplet splitting in this theory) 



then CP violation implies [lj that x must lie on one of the two complex solution 
branches x±(£), (£ G (—27.917, oo)). Since A, r\ are already counted as independent 
x+(£) counts for M/m. 

As in the case of the MSGUT one imposes the fine tuning condition DetTi = to 
keep a pair of Higgs doublets H(i), #(i) (defined via left and right null eigenstates of 
the mass matrix Ti) light [6J, [9J, [3], [1] . The composition of these null eigenstates in terms 
of the GUT scale doublets then specifies how much the different doublets contribute 
to the low energy EW scale symmetry breaking. In the Dirac mass matrices we can 
replace < hi >— > diV u1 < hi >— > diVd- The fermion Dirac masses may be read off 
the decomposition of 16 • 16 • (10 © 120 © 126) given in [7J El [3] and this yields(we 
have removed sin/3, cot (3 factors in h, g, f, r^ relative to our previous notation in the 
expressions for mass matrices [3J, 1] and written the result as h, g, f, fi 

V = {h + f + g) ; ri = — ; r 2 = — 

ai a 2 

f = (h-3f + (h-3)g) ; h ln/:>,f "' 



y d = (hh + hf + feg); r 6 



a>6 + iySctc, 
a 6 + iV3a 5 
a 6 + i\/3a 5 



y l = (hh - 3f 2 f + (f 5 - 3r 6 )<7); f 5 = , J (4) 

<7 = 2igJ-(a d + iV3a 5 ) ; h = 2\f2hoti ; / = -4W-«/a 2 

The Yukawa couplings of matter fields with 120 Higgs field give no contribution 
to the Majorana mass matrix of the superheavy neutrinos Va so it remains M AB = 
%\Pijabo ' ■ Thus the Type I contribution is obtained by eliminating Va 

W = -M D AB V A V B + V A m AB u B + -»• -M v A (I B ] u A v B + .... 

M^ 7) = -((mT(MTVO,4 B = 2«<&t£ (5) 

As shown in[TTll5] it is very likely that the Type II seesaw contribution is totally 
subdominant to the Type I seesaw. However the consistency of the assumption that 
it is negligible must be checked and quantified. Evaluating the tadpole that gives rise 
to the Type II seesaw (the 120 plet contributes new terms) one obtains pQ the Type 
II contribution to the light neutrino Majorana mass : 



2 

tI J = 16if A B < O- >= 16iV2f AB (i—=ai + iV6r]a 2 - VS(a 6 + iCa^a^-y-) 

\J I Mo 

= 2/c<2t£ (6) 



where Mo = 2(M + r](3a — p)). We have always considered Type I and Type II 
mechanisms together since in the NMSGUT one does not have any freedom to switch 
one off at will. 

As is apparent from the fermion mass formulae the coefficients «j, cti are critical 
for the phenomenology of these models. They are calculated by determining the null 
left and right eigenvectors of 7i (the [1, 2, ±1] mass matrix). The explicit (and rather 
complicated) expressions for the ai,at can be found in Appendix C of pQ. The only 
modification being that we work in a convention where a,\ = ol\ are both real so that 
the contributions of the 10— plet to all Dirac masses are real. 

From the formulae in[l] one can check that for real superpotential parameters and 
real values of x the «;, «j, i = 1...5 are real while «6, «6 are pure imaginary. Then it 
immediately follows[TJ that except the trivial 126 phase convention dependent phase 
values all other phases are zero and there is no CKM CP violation. The only way to get 
non trivial phases in the CKM matrix while retaining real superpotential parameters 
is thus for x to be complex. Thus we constrain the value of £ to be real and greater 
than -27.9 and choose the complex solution branch of the cubic equation which has 
a positive imaginary partpQ : this is an arc in the upper half complex x-plane with 
real value between about .75 and 2.25 while the imaginary part lies between about 
and l(see pQ for details about solutions of the cubic). 

3 Numerical search and anatomy of specimens found 

We use the 2-loop RG equations for the SM[30]and MSSM j3Tj[32] to extrapolate the 
central values of fermion yukawa data at the top mass mass scale (including neutrino 
mass splitting and mixing data), which were recently recalculated and summarized 
in [29l [26], past the MSSM threshold (taken to be 1 TeV and defined by expressing 
all soft masses in units of the gluino mass fixed at 1 TeV ) up to the MSSM 1 loop 
gauge coupling unification scale Mjj = 10 16 ' 3 GeV while ignoring, for the time being, 
complications such as righthanded neutrino thresholds. In the case of the neutrino 
masses we extrapolate the coefficient kab = M nu /(2v%) of the SU(2) x U(l) invariant 
dimension 5 operator [39], HO] that gives rise to neutrino mass when the Electroweak 
vev {y u = 174 sin (3GeV; for the SM sin/3 is 1) is substituted for the two electroweak 
doublets in it. 

At Mx we re-extract the input parameters i.e the 9 charged fermion yukawas, the 
4 CKM parameters, and the 5 neutrino mass data. Using the percentage uncertainties 
in fermion yukawas (for tan (3{Ms) = 50) at Mx = 10 ie ' 3 GeV recently re-estimated 
in [23], which include the large uncertainties due to large threshold effects 
on down type quarks and charged leptons at large tan (3 , we form a \ 2 func- 
tion comparing the central values with the values of the same parameters calculated 
using the formulae described in Section 2 and normalized by the errors mentioned. 
The 'down hill simplex' or "amoeba" algorithm of Nelder and Mead [25] is then used to 
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Parameter 


SOLN - 1 


SOLN - 2 


x 2 


5.71427 


6.79112 


/in 


8.95176 x 1(T 7 


2.69589 x 10~ 6 


^22 


0.000547776 


0.00053838 


^33 


0.196286 


0.19651 


/ll 


-2.06673 x 10- 7 


9.89865 x 10- 8 


/l2 


-2.37214 x 10~ 6 


2.02625 x 10~ 7 


/l3 


-6.34652 x 10~ 6 


5.56002 x 10~ 6 


/22 


-0.0000237101 


-5.28297 x 10~ 6 


/23 


-0.0000887439 


0.0000439064 


/33 


-0.000158592 


0.00011407 


#12 


-0.00592746 


-0.0258934 


9l3 


0.0392534 


0.0678888 


923 


0.428519 


1.81457 


X 


2.1294 


0.247415 


V 


-0.590395 


-0.619743 


P 


-0.0325812 


0.0682899 


k 


-0.000170967 


-0.000055478 


rho 


-0.0287934 


0.338748 


{m flt /(10 w GeV)} 


0.517858 


0.347328 


7 


-0.00131809 


-0.0003681 


7 


0.0421822 


0.04502453 


c 


-0.134801 


-0.057408642 


c 


0.291861 


-1.53557 


e 


0.848938 


0.596808 


{m BG /(10 16 )}Gey 


2.897 


.304376 


Ax 


1.07903 


1.60956 


"g 1 


14.7551 


1.06413 


Aa 3 (M z ) 


-0.00132645 


-0.015022 


{M^7(10 12 )}Ge\/ 


{39.5853, 3.96846, 0.00122388} 


{48.476,7.463,0.001046} 


{M / l/ / /(10- 12 )}e\/ 


{25.08,2.51,0.000775} 


{0.0392, 0.0061, 8.48 x 10~ 7 } 


M^ rc (Gey) 


{82.002,6.485,0.01606} 


{82.257,6.468,0.0164} 


M„{meV) 


33.934,6.146,1.819 


33.748,5.96,0.9995 


{Eigenvalues[f}} / (10~ 5 ) 


-20.283, 2.033, -.00062 


12.8704, -1.9815, -0.0002779 


M H (units of m) 


-0.00123516 + 0.000481764i 


0.000633884 - 0.000318149i 



Table 1: Two sample sets of NMSGUT couplings which pass the criteria of Section 2. 
Note how the strong hierarchy of M v " weakens the strong hierarchy of (M£, rc ) 2 to re- 
produce the weak observed light neutrino hierarchy via the Type I seesaw mechanism. 
Note also the utter irrelevance of Type II neutrino masses. 



troll the 23 dimensional hyperspace of NMSGUT parameters searching for local mini- 
mae of the x 2 function. While searching we imposed penalties to keep all the Yukawas 
and GUT superpotential parameters within the perturbative range. The programs 
were written in FORTRAN and run on PC's for 10 6 iterations for each random start- 
ing point (or less if it was clear that the algorithm was stuck in a local minimum with 
large x 2 '■ a ^ which point the program was restarted). Typically our algorithm found 
candidate' solutions' every 20 program iterations or so with 4.0 < x 2 < 8 (but never 
any lower) . 

The striking features of the candidate solutions (besides the inability to pass below 
X 2 = 4.0) were : 

• Neutrino masses and large mixing angles were accurately fitted. 



• The values of the 126 couplings were so tiny that the three righthanded neutrino 
masses were much lighter than Mx and spanned the range ~ 10 8 — 10 13 GeV, 
just as conjectured by us earlier [3j [121 03] • Note that this is just the range 
required by most (Type I) Leptogenesis models. 

• The bulk(over 85%) of the residual x 2 typically consisted of errors in the d,s 
quark yukawa couplings with a further 10% or so coming from the bottom mass. 

• The overall scale parameter of the MSGUT symmetry breaking which suppresses 
the neutrino masses was typically in the range 10 15 — 10 16 GeV so that the role 
of the tiny 126 couplings in generating viable neutrino masses is manifest. 

We then subjected the candidate solutions to a highly demanding ('ouroborotic') 
consistency test between the neutrino mass fitting specified scale parameter m (since 
this fixes the scale of the right handed M p mass matrices and thus the light neutrino 
masses) and the value of m determined [1R EJ [1] by the RG flow via the specification 
of Mx- We required that they should lie within a factor of 5 of each other. In view 
of the many approximations and short cuts adopted in this initial survey of the huge 
NMSGUT parameter space this requirement is quite stringent but we even found 
solutions which agree to within 15% for these crucial numbers. We also demand that 
the correction to as(Mz) due to the GUT scale threshold correction effects evaluated 
using the complete NMSGUT superheavy spectra[9j [lj [36] should lie[35] within the 
range —0.017 < Aas(Mz) < —.004 and that the unification scale correction should 
raise [I] the unification scale (scale of d = 6 baryon violation). We have relaxed the 
upper limit to —0.001 for the case of Soln-1 since it was the first 'good' solution we 
found: subsequently we always use —0.004. These three demanding requirements ( 
which can perhaps be thought of as adding 3 equations to the still deficient set of 18) 
killed practically 95% of the randomly generated solutions immediately. Nevertheless 
among every 15-20 candidates we would find one that seemed to be suffering only 
from the 'ds-br disease' : which we propose is rather a legerity- adiposity attained 
when crossing the supersymmetric portal. 

10 



We then extrapolate these 'rough nuggets' back down to the summarily assumed 
Susy breaking scale Ms = 1 TeV. Only the d, s quarks were seriously in error while 
b, t require minor revisions upwards from their standard model values. We use the 
freedom to choose Susy soft parameters and the elementary formulae matching the 
threshold corrected MSSM down/charged lepton yukawas to their SM analogues [23] 
to correct the discrepancies. These corrections all have [23] the form(for the i'th down 
type fermions only): 



y? UT (M s ) 



.SMi 



{cos /3 + ei(mp.,mjj j, M!,M 2 ,M 3 ,ij,, A t , ai,a 2 ,a 3 ) sin/?} 1 (7) 



yr J (M s ) 

The explicit forms of these formulae and the consequent larger uncertainties at 
Mx were summarized recently in[23] and we used those uncertainties while locating 
candidate parameter sets. The equations given by (23] following (21] for ej are 

2as /i 



ef + ef + ef + e y S ib , with 



G 



B 



w 



WW* Ha ( u Q<> u v > 



i 



167T 2 
1 



'3 

,/2 



9_M] 
6 fx 

3g 2 M; 



1 -(H 2 (v^,x 1 ) + 2H 2 (v^ 



Xi 



/2 

9 V_ 

9 M x 



H ^QV W di 



16tt 2 2 n 



H 2 {v .x 2 ) 
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if = m 2 -/ Ml, Vf = rrif/fr, wj = myM 2 , x\ = M\j ' u 2 and x 2 = M$//jr for 
i = d,s, b. All mass parameters are assumed to be real. The correction e y to d, s- 
quarks are neglected because of the strong hierarchy of the quark Yukawa couplings. 
H 2 is defined as 
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H 2 is negative for positive x and y and l-f^l is maximal, if its arguments are minimal, 
and vice versa. It is a rather smooth and bounded function with typical values 
~ 1(T 2 - KT 1 . 

The last term in equation f[T3|) causes an enhancement by a factor of —9 compared 
to the corresponding term in the quark sector in equation ([9]). The r-lepton Yukawa 
coupling does not have a relevant correction of the e y type because the corresponding 
vertex correction is suppressed by the heavy mass scale of the right handed neutrinos. 
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for i = e, fi, t. 

Taking the sign of M 3 positive by convention we see that the dominant as de- 
pendent corrections to the down and strange yukawas, which must be lowered at Ms 
from their SM values since the fit values at Mx are so small, will be positive only if 
/i is positive and large relative to the Gluino mass. Conversely the bottom yukawa 
coupling in the MSSM must be raised at Ms from the standard model value and 
hence the value of A t should be negative and appreciable relative to fi (note that 
while y t is ~ 1 there is a relative factor of 6 overall in the formulae for eQ,e y ). On 
the other hand, for positive values of Mi the terms in lepton q ' are both positive 
and negative. Thus various solutions with both signs for M\^ are possible. Since 
the interpretation of Gaugino phases is fraught with complexities we usually prefer 
to display only solutions with positive signs for both M 12 . 

However the final fit is at Ms- The 19 (!) free parameters in these formulae 
make the task of trimmming and padding the the down type fermion Yukawas to 
match at Ms facile using an appropriately tuned downhill simplex program. We 
obtain fits accurate to a few parts in 10 16 for the d, s, b, e, /i, r fermions. This clearly 
indicates that the sfermion parameters are underdetermined and there is, in principle, 
ample room to accommodate constraints from Flavour violation, Baryon Decay and 
Unification. 

An interesting point in this regard was that we did not prejudge the value of 
tan/3 but used it as another fitting parameter. Since only the low energy values 
of the fermion yukawas and Weinberg neutrino mass operator coupling k have any 
sanctity while the consistent GUT data sets are taken at face value perhaps our proce- 
dure is not completely wild. A further consistency check is then that the T 3 ^ = +1/2 
quark masses and the neutrino masses should also lie within experimental errors and 
we retain only such soft parameter sets for which these errors are within the stan- 
dard experimental errors evaluated at m t op{m t0 p). Although we also find solutions of 
the fit at Ms which give extremely accurate fits to the lepton masses comparable to 
the extreme precision with which they are known at Mz it should be kept in mind 
always that our results carry systematic errors due to the omission of various ef- 
fects such as right handed neutrino thresholds (important due to their light masses), 
due to the mass splitting among Susy partners at M s and due to non-diagonality 
in the supersymmetric threshold corrections. The right handed neutrino thresholds 
will particularly affect the neutrino parameters and so they should be allowed some- 
what larger margins of error. The first of these is quite amenable to correction and 
the second can at least be estimated and the soft spectra found subjected to self- 
consistency The last mentioned theoretical results can be easily calculated from the 
same diagrams^ by inserting the appropriate flavour structure and diagonalization 
and this is now urgently required for the next stage of our program. However the 
results depend on a large number of unknown sfermion mixing parameters and it 



1 C.S.A is obliged to Ts. Enkhbat for a discussion on this subject. 
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will be necessary to make simplifying assumptions in order to proceed. Luckily, the 
constraints on Flavour violation will strongly suppress most of the mixing angles, the 
problem is mainly to disentangle the various dependencies. One may also extend the 
number of targets of the random fit at Ms by optimizing also the up quark masses 
and the neutrino mass squared splittings : however the only parameter that affects 
these is tan/3 and we found it opportune not to overburden this formerly quite free 
number with too many sartorial duties ! When one includes non-diagonal threshold 
corrections(which will presumably affect the CKM matrix as well as lead to flavour 
changing processes) the entire sfermion sector including the the sup and scharm quark 
masses will enter the calculation [49|. 

In general, due to the the structure of the threshold correction formulae and the 
opposing requirements of strong corrections of charge -1/3 quarks and minor fine 
tuning of charged lepton masses, the fits are easier if one allows oneself the somewhat 
moot luxury of negative gaugino mass parameters. Since the role of phases of the 
gaugino masses is obscure we will allow both signs but only real masses. When 
searching for the soft parameters which would cure the d-s-b quark mass disease we 
imposed various commonly accepted requirements (to within 10% or so) on the soft 
parameters such as the mass ratios M\ : M2 : M3 :: 1:2:6, degeneracy of 
the first two generation sfermion masses, and the constraints from RG evolution on 
the ratios of the masses of the sfermions of the first two generations to the gaugino 
masses [27] : 

\fh TtL /M 1 \ > 0.9 ; \rhg, Q /M 3 \ > 0.75 (15) 

However the conventional wisdom that the sfermions of the third generation will 
be much lighter than the first two generations of squarks due to strong renormaliza- 
tion by top yukawa driven effects seems directly challenged by our spectra which seem 
always to include a large mass for at least one of the third generation sfermion mul- 
tiplets at M s . Since one can obtain solutions where this could be any of f, b, i, Q 3 , L 3 
it may still be possible to arrange that this happens in a way consistent with the 
expectations from the RG flow : for instance if the sfermion with the largest gauge 
contributions to its mass corrections(e.g (^remained massive while the other third 
generation fermions are driven to lighter values as usual. If however we demanded 
that all the sfermions of the third generation are lighter than the corresponding first 
two generation fermions then we have, so far, failed entirely to find any fits at all ! 
This feature thus marks out such sfermion spectra as linked to the inner logic gov- 
erning the fermion mass pattern in the NMSGUT and we proffer it as a prediction 
regarding supersymmetric physics that may be tested by LHC or its successors. The 
trilinear soft parameter A t for the third generation and the /1 parameter from the 
superpotential play critical roles in the simultaneous correction of the the b mass 
downwards(as is generally required since the b quark yukawa tends to drift towards 
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the r yukawa in SO (10) fits rather than stay at the - upto 30%- lower value indicated 
by the MSSM RG flow at large tan/3. Although we have set the gluino mass M 3 mass 
to 1 TeV for definiteness since the threshold correction factors depend only on ratios 
of masses it is clear that one may need to consider a range of values in general. 

Given the large value of u one may wonder how EWRSB (always a bit difficult 
to implement in large tan/3 models ) will work here. This depends strongly on the 
assumption for the Higgs soft masses and since it is quite possible to envisage scenarios 
where these masses are tachyonic there is no reason to be overly pessimistic in this 
regard. The drastic constraint of suppression of d — 5 mediated proton decay is 
a further dragon that the parameter sets found will need to face. Nevertheless the 
very fact that they came so far while satisfying all structural requirements is, in our 
view, highly remarkable in such a constrained and structurally unique theory as the 
NMSGUT. 

The parameter sets found so far display some characteristic preferences which 
mark them out as possible signals of NMSGUT driven behaviour at LHC scales. In 
Tables 1,2,3,4 we give example NMSGUT parameter sets which satisfy our criteria. 
However this is done only as a proof of principle and feasibility and the 
typical form of non pathological exemplars of our fitting program. It 
would be highly premature to subscribe to any one parameter set as very 
specifically demanded by NMSGUT fits. Although the fits truly involve only 
the yukawa couplings and Weinberg operator coefficient kab (since the low energy vev 
is supplied by hand on both sides of the mass formulae) we display the values of the 
yukawas times the two loop run up of the MSSM Higgs vevs just to make the numbers 
somewhat recognizable via their recognizable origin the low energy mass values. We 
defer the study of the properties of the sfermion spectra found at the GUT scale vis 
a vis the requirements of unification to a later more extensive work[49]. Clearly the 
highly structured solution will lead to differences from the standard SUGRY scenarios 
and require generation dependent soft terms at Myp], 

3.1 Dissecting Specimens 

The basic difficulty in building a grand unified theory of fermion masses is the 
achievement of sufficiently large neutrino masses and the reconciliation between small 
CKM mixing angles and large PMNS ones. So it is interesting to dissect the solu- 
tions we have found to see how exactly they circumvent the difficulties typically 
encountered [TTJ [3j [121 02] • The roles played in the effective matter fermion Yukawa 
couplings Y u,d,e , kab by the Yukawa couplings h,AB,gAB, Jab of the three FM Higgs 
representations are, so to speak, the vital organic functions of the NMSGUT and 
their dissection reveals the integral unity and 'intelligent design' within the fermion 
mass pattern. By killing the Yukawa couplings individually or in groups and noting 
the resultant effect one easily uncovers their roles. 
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Parameter 


Target 


Uncert. 


SOLN1 


SOLN2 


m t {GeV) 


76.6 


3.06 


75.69 


75.96 


m c (MeV) 


194.97 


32.2 


207.31 


193.477 


m u 


.4 


.15 


.41 


0.39 


m b 


1186.8 


616.0 


1621.31 


1626.01 


m s 


17.82 


8.4 


4.47 


4.92 


m d 


.931 


.54 


.24 


0.25 


m T 


1493.37 


283.74 


1693.83 


1698.87 


"V 


74.989 


11.248 


75.01 


76.65 


m e 


0.355 


0.0533 


.35 


0.35 


sin 6\ 2 


0.22098 


0.0016 


0.221004 


0.220974 


sin 6*23 


0.0349 


.0013 


0.0347803 


0.0349877 


sin 6*13 


.00297 


.0005 


0.00294978 


0.00305354 


Scp 


60.026 


14.0 


52.0769 


64.53 


sm 2 9[ 2 


.30029 


.06 


0.298219 


0.288079 


sin 2 9^ 


.465 


.14 


0.467765 


0.326723 


sin 2 #f 3 


.0053 


.03 


.0108116 


0.0292243 


m 2 2 (xlO~ 5 eV 2 ) 


3.46 


.37 


3.44621 


3.4533 


m 2 23 (xlO~ 3 eV 2 ) 


1.11 


0.22 


1.11379 


1.10343 


x 2 





18.0 


5.71427 


6.79112 


xl 





1 


1.58 


1.536 


xl 





1 


2.51 


2.351 


xl 





1 


0.49 


0.508 


Xdown 





1 


4.58 


4.395 



Table 2: Values of the fermion data at Mx achieved by the sample coupling sets 
of Table I. The values v(M x ) = 136.729 GeV, tan (3(M X ) = 46.819 obtained by a 
complete two loop extrapolation to Mx of central values at Ms [29j have been used 
to write the notional masses at the GUT scale instead of the Yukawa couplings and 
neutrino mass generating Weinberg dimension = — 1 couplings hab- 
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Parameter 


SOLNl - Set! 


S2 


S3 


SOLN2 - SI 


S2 


S3 


Mi 


0.17876 


0.1687 


-.1786 


0.165479 


0.1759 


0.1677 


M 2 


0.3260 


0.3315 


-.3703 


0.33562 


0.3326 


0.3299 


M 3 


1.0000 


1.0000 


1.0000 


1.0000 


1.0000 


1.0000 


A t 


-6.7908 


-27.1390 


-18.9522 


-1.41918 


-7.4401 


-3.4534 


n 


22.6655 


17.8729 


16.4992 


9.65223 


14.4994 


21.4061 


™<)i 


3.7317 


2.2145 


1.16261 


0.82551 


0.8975 


2.7146 


m Q2 


3.4621 


2.1701 


1.2201 


0.90581 


0.8949 


2.9140 


™Qz 


5.2620 


18.0030 


5.209 


1.125382 


20.0218 


1.1453 


rh Ll 


0.8036 


0.6337 


2.0285 


0.812268 


0.6688 


1.3384 


m L2 


0.7623 


0.6406 


2.0974 


0.897360 


0.7391 


1.4735 


m L[i 


0.3000 


0.4535 


1.2242 


0.3634 


0.4303 


3.6479 


m d 


0.8969 


1.6780 


2.6808 


1.85571 


2.6043 


1.7563 


m s 


0.9840 


1.5881 


2.4728 


1.80036 


2.6405 


1.6093 


m b 


27.6858 


1.7006 


7.2372 


29.29089 


16.9781 


29.9958 


m t 


0.7858 


0.3673 


1.5447 


0.421398 


0.3259 


0.5962 


m e 


1.4492 


1.3102 


4.6487 


0.39730 


0.6527 


0.5595 


«V 


1.5055 


1.3102 


4.3451 


0.436624 


0.6796 


0.6076 


fh T 


7.3105 


6.1334 


0.6736 


2.02467 


1.6692 


0.5264 


tan(/3) 


50.3352 


50.0910 


45.7251 


51.1654 


54.3213 


52.0464 



Table 3: Values (TeV) in of the soft susy parameters found by the threshold fitting 
at scale Ms- Note the large values, relative to the isoquantal first two generation 
fermions, of at least one of the third generation sfermions in each column. 



16 



Parameter 


SM(M<?) 


Uncert.(%) 


Ach.value. 


Ach.value 


Ach.value. 


Ach.value 






(SM) 


SolA(Sl) 


SolA(S2) 


SollO(Sl) 


SollO(S2) 


m b 


2394.77 


2.1 


2394.77 


2394.77 


2394.77 


2394.77 


m s 


45.9467 


28.3 


45.94676 


45.9466 


45.9466 


45.9466 


m d 


2.43876 


42.2 


2.43876 


2.43876 


2.43876 


2.43876 


m t 


145873.0 


1.65 


147472.0 


147472.0 


147303.0 


147306.0 


m c 


520.549 


13.3 


526.254 


526.253 


525.652 


525.664 


m u 


1.07639 


36.5 


1.08819 


1.08819 


1.08694 


1.08697 


m T 


1734.45 


1.12 x 10~ 4 


1734.45 


1734.45 


1734.45 


1734.45 


"V 


102.014 


8.88 x 10~ 6 


102.014 


102.014 


102.014 


102.014 


m e 


0.483143 


8.6 x 10~ 6 


0.483143 


0.483143 


0.483143 


0.483143 


m 2 2 (meV 2 ) 


81.3871 


3.75 


83.1851 


83.1845 


84.5943 


84.6016 


m2 3 {meV 2 ) 


2543.34 


8.0 


2630.3 


2630.28 


2510.63 


2510.85 


sin 9\ 2 


0.2210 


.7 


0.221019 


0.221019 


.220993 


.220993 


sin #13 


0.00345 


3.7 


0.00351868 


0.00351868 


0.00364475 


0.00364475 


sin 0f 3 


0.04059 


15.5 


0.041483 


0.041483 


0.0417563 


0.0417563 


5cp 


60.024 


23.3 


52.08 


52.08 


64.535 


64.535 


sm 2 6[ 2 


0.31 


7.4 


0.314852 


0.314852 


0.2904 


.2904 


sin #{3 


0.005 


<0.03 


0.012641 


0.012641 


0.034 


0.034 


sin 2 6 23 


.500 


4.0 


0.5082494 


0.5082494 


0.3637 


0.3637 



Table 4: Values of the precision data at m t {m t ) extrapolated to Ms = 1 TeV compared 
with the values achieved by the NMSGUT parameter sets listed in Table 1 extrapo- 
lated down to Ms and translated to Standard Model values by applying threshold cor- 
rections using Susy spectra and tan/3 selected using a downhill simplex optimization. 
Quark masses are in units of MeV. The mixing values obtained from the runddown 
from Mx have been listed unchanged since only diagonal threshold corrections were 
applied. Systematic theoretical errors of the GUT parameter estimates have not been 
included in the error estimates which can be only notionally as precise as quoted for 
the leptons. 
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The most dominant coupling (though not, generally, the largest) is h^ which di- 
rectly and drastically affects every mass and mixing angle and may justifiably be 
called the lynch pin of the yukawa coupling pattern. The third generation masses can 
be fit to within 5% by this coupling alone : this is precisely the content of '6 — r — t 
unification' in this context : even though the Yukawa couplings Y u ^ e (M x ) themselves 
may differ by as much as 30%. The next most important( and largest) coupling is 
#23 whose value is generally larger than /133 perhaps to compensate for the fact that 
it enters the mass matrix non-diagonally. It is crucial for the #23 mixing and thus 
also for second generation masses: specially the u mass. The importance of the cou- 
plings Jab arises from their smallness : since they control the mass and couplings 
of the superheavy right handed neutrinos which leave behind a characteristic signa- 
ture in the Type I seesaw masses (where they enter in the denominator) the limit 
Jab — > is singular i.e decreasing Jab changes the effective theory (MSSM with 
massive left handed neutrinos) drastically. The couplings hn, /i 2 2, #12, 9i3 are relevant 
for the masses and mixing of the first and second generations. 

If we first set j ab — then clearly all neutrino masses and mixings become 
unphysical. However one finds that the fractional changes in the charged fermion 
masses and mixing are changed by less than one part in 10 4 — 10 6 . Thus the 126 
coupling is utterly irrelevant to the charged fermion sector in our fits. 
This feature clearly marks out the solutions we have found as totally 
distinct from the class of generic fits explored in [T1^[T5]. Moreover our Jab 
eigenvalues (and thus the right handed neutrino masses) are strongly hierarchical 
(10 2 : 10 1 : 1CT 3 ); so much so that they are effective in compensating the large 
hierarchy in the neutrino Dirac masses(10 2 : 10 : 10~ 2 ) (which is amplified by the 
square of the neutrino Dirac mass in the Type I formula) to produce the rather 
weak (15 : 3 : 1) hierarchy of the light neutrino masses. This is again crucially 
different from the generic SO(10) solutions of [HI US] which have a hierarchy of 
at most 300 : 60 : 1 for the right neutrino masses ( which are also much larger). 
There [HI [15] the dilution of the Dirac mass hierarchy relies on a strong boost of 
the coefficient r D of qab in the Dirac neutrino mass matrix M^ TC to huge values 
~ 10 4 — 10 5 . Our parameters simply never venture towards such large rr> '■ they 
always produce values of the coefficient ro below 500. In fact even if one removes 
constraints that the parameters remain perturbative it proves extremely difficult to 
force the downhill simplex 'amoeba' -in its unicellular wisdom !— to crawl towards 
the high r^ region. Of course we have already shown [2] that generic numeric fits 
carry no direct implication for the NMSGUT since their parameters are constrained, 



but still the generic analysis propagated a scare [18] that the small 126 scenario was 
utterly flawed and then compounded it by the confusing conclusion that small (but 
not tiny) Jab solutions were in fact capable of finding accurate fits. In any case this 
history underlines once again [TR |3] the emphatic non sequituur with which explicit 
UV theories must confront such generic no-go claims unless they are founded on 



something more than numeric projections. 

In JT2J H3] we argued that the 2-3 sector should be regarded as the core of the 
fermion mass hierarchy and that the complete hierarchy could and should be under- 
stood semi-analytically by perturbing around the couplings ^33, #23 i-e by taking these 
couplings as 0(1) and the others as order 0(e n ),e ~ 0.1, n > 2. The solutions we 
have found by random search bear this out fully. The couplings in the 12 sector are 
hierarchical and suppressed relative to the 23 sector : h 2 2 ~ ^h 33l h\\ ~ £ & h 33) g\ 3 ~ 
e 923 , 9i2 ~ e2 fi , 23- If we keep only the couplings h 33 , g 23 and set the rest to zero one finds 
that Tn tt b,T,fjLi $23 are s tih within 1.0% or so of their final values. Restoring h 22 ~ e 3 
restores also m CjS to within 5% or less of their final values and improves #23, but the 
CKM parameters apart from 6 23 are still 100% or more different. The switched-on set 
{^33> <7i3, ^22} gives good values of t, c, b, s, r, 13 but not the rest. Similar considera- 
tions apply to {h 33 ,g u , h 22 }. Finally restoring #12, #13 restores the rest of the CKM 
parameters i.e 6\ 2 , 0f 3 , 5cp to within .05%. Thus the qab are directly correlated with 
the corresponding 6ab angles. We can now pursue [50] our earlier proposal [12, [13] 
to understand the fermion parameters by an expansion around the 23 sector in a 
single hierarchy parameter with confidence that the solutions exist and have the re- 
quired form. The phases and coefficient values we have found can be checked for 
their compatibility with such an expansion and thus a qualitative understanding of 
the hierarchy arrived at. 

Of course we must keep in mind that the values of the Yukawa couplings we find 
match the SM precision couplings extrapolated to Ms only after applying the large 
corrections at Ms due to large tan /3 that are the underpinning of our work. Thus the 
fermion masses explained above are themselves significantly different from the ones 
that we earlier[12j [13] had only partial success in describing via the e— expansion. 
Indeed it precisely this re-scaling downwards of the target ma, m s in the fit at Mx 
downwards in the hierarchy by a factor of e that we anticipate will allow the determi- 
nation of excellent fits to the large angles of the lepton sector. Analytic (perturbative) 
confirmation of this proposal is thus an important topic for future work[50j since it 
will yield a more robust understanding of the innards of the hierarchy. 

4 Nucleon decay 

As was already shown generically inpQ and can also be seen from Table 1, viable 
RG flows in the NMSGUT have Mx raised by up to two orders of magnitude. Thus 
d = 6(i.e gauge boson mediated) proton decay will be suppressed by up to 8 orders of 
magnitude relative to the already long life times ~ 10 36 yrs corresponding to the one- 
loop unification scale 10 16 ' 3 GeV. Thus question of the rate of d = 5 operator mediated 
baryon decay in the theory at hand is clearly the critical one. Since the decay rates 
in these channels will depend sensitively on the soft scalar spectrum and mixing they 
could provide a welcome additional filtration of the parameter sets that pass the 
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other criteria. Of course since the threshold analysis so far is all performed in the 
diagonal approximation the mixing parameters of sfermions of the same SU(3) x U(l) 
quantum numbers are arbitrary and one could at most derive upper bounds on these 
mixings once the variation possible due to sfermion spectra had been mapped. For 
the sfermion masses on the other hand we have already used the freedom to choose 
all but the right handed u, c squark masses. Of course -given the non-uniqueness 
of the parameter sets we obtained- it is possible that a more general analysis could 
later impose the requirements of fixing the down fermion masses, of Unification(of 
soft masses), obtaining acceptable proton decay partial widths and flavour violation 
due to loop effects simultaneously. 

For the moment however we will simply treat the u,c squark parameters, the 
squark mixing, soft trilinear couplings^^ ), Higgs scalar B parameter and the 
Higgs soft masses (-C S oft = BHH + m 2 H \H\ 2 + m 2 j\H\ 2 + ...) as the only remaining ad- 
justable parameters in our search for completely consistent data sets in the NMSGUT. 
Actually, at large tan/? one can quite plausibly fix the last three parameters via the 
Electroweak symmetry breaking conditions ; in any case they affect sfermion mass- 
mixing structure and thus nucleon decay only through their contributions to the the 
two loop beta functions. The dominant vev v u may be fixed (to 0(cot/3)) from the 
leading approximation to the scalar potential which retains only v u thus fixing one 
soft parameter out of B,m 2 1 ,m'g. The minimization conditions will eliminate two 
more. 

The suspicion may be that like others jlSJ 46J this GUT too will be hard put to 
account for evasion of the quite stringent experimental bounds on the commonly ex- 
pected modes in supersymmetric GUTs. A proper analysis of this question requires 
extensive calculation of the large arrays of Wilson coefficients: which we have cur- 
rently undertaken. However we can get a preview of the issues, techniques and results 
by considering the channel N — > K + u T since they are the dominant ones at large tan/3 
anyway [421 [28]. 

In [9], [1] we calculated the effective d — 5 operators for nucleon decay in terms 
of the couplings h^B-, Iab^Qab and the (inverses of) the mass matrices for the three 
different triplet types relevant to proton decay in this model. Not only the i[3, 1, — |]ffi 
t[3, 1, |] triplets (square brackets enclose MSSM gauge quantum numbers and the 
subscripts refer to the multiple copies arising from the various grand unified Higgs 
present in the theory : whose spectra have been completely evaluated [7J [H El El 
[TUl [3l []]) but also in the novel ones from the P[3, 3, ±|] and K[3, 1, ±|] multiplet 
types which can contribute to baryon violationpQ H7j. In the case of the 126 the 
Pi, K\ C 126 multiplets couple to matter fermions but P\,Ki C 126 did not. The 
120 however contains P2, Pj and K 2 ,K 2 . Since these mix with Pi, Pi and Ki,Ki, 
a number of fresh contributions to Baryon decay appear(although they are usually 
suppressed relative to the usual triplets by their larger masses). 

On integrating out the heavy triplet Higgs supermultiplets one obtains [9l [1] the 
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following additional effective d — 4 Superpotential for Baryon Number violating pro- 
cesses in the NMSGUT to leading order in mw/M X - 

Wf}f^° = -L ABCD (-€Q A Q B Q C Ld) ~ RABCD{^AU B U C d D ) (16) 

where the coefficients are 



"ABCD — S\ hABhcD + <$l h-AsfcD + $2 f AB^CD + S 2 JabJcD 

— <S 1 h AB gcD — S 2 Jabc/cd + v2(V~ ) 2 (jacIbd 

- (V- l ) 2 2 ~g A c~g B D (17) 



and 



Rabcd — S\ h A BhcD — S\ h AB fcD — S 2 Jab^cd + S 2 JabJcd 

- iV2S A l f AB h CD + iV2S A 2 f AB f CD 

+ ^6 QABhcD — iS 7 g AB h C D — <S 6 cjabIcd + ^7 9abIcd 

+ iS^~h AB ~g C D ~ iS 2 7 f AB g C D + V2Sjf AB g C D 

+ iS e 7 g AB g C D + S 7 7 g AB g C D - V2(fC' l ) 2 f AD g BC 

- {K,- x ) 2 2 ~g AD ~g BC (18) 

here S = T _1 and T is the mass matrix for [3, 1, ±2/3]-sector triplets : W = 
tT/tj + ..., while 

h AB = 2V2h AB f AB = 4y/2f AB g AB = Ag AB (19) 

The carets on the fields denote that they are in the Hab diagonal basis(at Mx) in 
which we first determined our parameter sets by a random search. The MSSM fields 
do not carry carets and are unitarily related to the GUT fields [TTj: 

Qab = Q*a'aQa' ', La = C* A , A L A i ; 

&a = (y u y A , A u A , ■ 2 A = (v d y A , A d A , ; i A = v e \, A e A , (20) 

The fields in terms of which we actually require the decay operators are those of 
the MSSM y u ' e -diagonal Y d diagonal-modulo-CKM basis. So one passes to that 
basis by the above standard bi-unitary transformations but with all the matrices 
determined (upto conventional phases) by diagonalization of the caret basis. The 
determination of the the right handed mixing at Mx underscores the completeness 
of our parameter sets. The values of these coefficients in the MSSM basis are thus 
easily computable and we spare the reader the explicit sample expressions . 
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A proper analysis, now under way|49j. requires one to integrate the renormal- 
ization group flow of the above set of coefficients down to the scale Ms along with 
the rest of the MSSM hard and soft parameter flow and there again pass to the 
MSSM Y u,e,d diagonal basis appropriate to that scale by rediagonalization. One 
will then dress the 2-fermion-2-sfermion operator contained in the supersymmetric 
AB ^ 0, d = 4 effective superpotential by gaugino/Higgsino exchange to obtain the 
effective 4 fermion operators for Baryon decay at Ms and finally renormalize these 
via SM RG flow to the nucleon mass scale. Although important for precise numbers 
these re-scalings do not add anything new as long as one has already allowed for an 
arbitrary scalar mixing matrix at Ms- Thus we can estimate the RRRR channnel 
decay rates for the first solution set simply using the {Rabcd(Mx)} array values(in 
the MSSM basis(without carets) at M x ) given explicitly above but for the rest us- 
ing the standard [HI H2J [28] chiral Lagrangian based formulae for converting from 
the quark to hadron states. The sfermion fields of the d — 5 effective operator in 
the Lagrangian are in the MSSM Y u,d ,e diagonal CKM basis, which is not the one 
that diagonalizes the sfermion mass squared matrices (which are in fact 6 x 6 if even 
0(M w /M s ) effects- which mix SU(2) L singlets and doublets -are retained). Thus 
their contractions (propagators) introduce the unitary diagonalization matrices (Wj 
of the the sfermion mass squared terms when rewriting these dressing contractions 
of the 2-scalar line, one fermion line dressing triangles in terms of the sfermion mass 
eigenstate propagators. 

The one loop triangle diagrams that dress the d — 5 chiral vertex from the GUT 
can be written in terms of a standard function : 

f(M, mu m 2 ) = — o o \,J n-r± 9 2 ^J n— |) 21 

JK ' ' ; m\-m\ K m\-M 2 M 2 m\- M 2 M 2 ' y ' 

Thus to account for the complete RRRR amplitude one need compute only the 
folding of this function - for a given sfermion spectrum- with the Rabcd coefficients 
estimated as approximately those given above at Mx and the unitary sfermion di- 
agonalization matrices. For the amplitude for the dominant p -^ K + v T channel this 
involves 

A = a{R ext {l + ^±*?-) + Rj-^-) (22) 

on Stub 



where 



R ext = MR* 3131 - ^311)^33^22 J2f(»,K,™l B )w™w? B 



B=l 



(23) 



R mt = ^((^ 132 - ^12)^3^1 E /(^ ™'. O^b^) (24) 
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U B J 

B=l 



The definitions of the baryon mass m# and the proton matrix elements a p , D, F, B 
can be found in [42] and we have indicated only the dominant 3 generation contribu- 
tions for the lepton side of the triangle diagram. 

Since all ingredients are in hand one obtains with Ai = 1.43; F = 0.44; D = 
0.81; a p = — .003GeV 3 , with a standard CKM type parametrization for the mixing 
matrix W u in terms of mixing angles (we assume the phase is zero for simplicity) 
(9 = 0i2,0 = $13, X = $23) , and with the free squark masses taken to be m s = 
.98mj, rhc = .99mf, for the proton life time in years and all masses in MeV 



_ : , w -„„ ( 32irm;/j) (3F + D) 2m p D 

ml-m\y MRext{l + mp 3m B )+Rmt j23^j l 



r p = 2.lxlO- 29 yr-MeV±- p -^-\a p {R ext {l + m p ^- >-) + R v 



p 

1I8 „,„„ xx /„„ n 2 



= 5.02565 x 10 yrs x (cos </> cos x sin x 

+ 1. 00061 (cos x sin 9 + cos 9 sin sin x) (cos cos \ sin — sin 9 sin x) 

— 1.00031 (cos 9 sin % + cos x sin sin #) (cos # cos x — sin sin 9 sin x))~ 2 

= 1.5 x 10 3 Vs (25) 

where the last equality is obtained at the angle values {0 = 0.001, 9 = 0.2, x = 0.002}. 
With only moderately different values of the parameters one can increase or (more 
easily) decrease this life time by factors of up to 10 5 or more. Thus we see that 
the result of the answer is extremely sensitively dependent on not only the sfermion 
masses and their degree of degeneracy but also on their mixing. The mixing of the 
sfermions will be limited also by their effect on the effective CKM mixing as well as by 
all the constraints coming from an analysis of supersymmetric quark Flavour viola- 
tion through box diagrams [43 j as well as of supersymmetric lepton Flavour violation 
though [44J. In our approach the soft susy parameters are thus doubly constrained : 
by the constraints on flavour changing processes and by accurate low energy fit to 
fermion masses using GUT parameter sets compatible with the high energy extrapo- 
lated data but with enlarged errors. Moreover we must still incorporate the effects of 
right handed neutrino thresholds in the RG flow of the candidate parameter sets down 
to Ms- An in extenso calculation of all Baryon decay widths in this model using an 
assumed soft susy parameter spectrum as used in e.g [42] would simply have no value 
in the current context and it would be inane to pursue or demand it. Thus a proper 
analysis will be lengthy, indeed it is an extensive research program that requires many 
innovations to be able to handle so many diverse constraints on parameters involved 
in fits to low energy data through highly non linear relations comprising not only 
highly non linear definitions of the effective MMSM Yukawa couplings (via the fine 
tuning of the doublet mass matrix) but also the RG flow between Mx and low en- 
ergies and the experimentally bounded but otherwise completely unknown sfermion 
data as well. This requires an order of magnitude increase in the scale of integration 
of our algorithms to enable searches that are automated from beginning to end. We 
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are still evaluating if this is worth attempting at this stage. We therefore leave a 
detailed study of the possibilities for a sequel[49j. 

5 Discussion 

In this paper (which is the second in a triplet dedicated to formulating the falsifying 
constraints to which the NMSGUT can be subjected on the basis of its structurepQ, 
fermion mass fits HE] and exotic processes [19]) we have indicated a feasible route 
to pursue if the New or Full Minimal Supersymmetric SO(10) GUT is to cross the 
final hurdles and become vulnerable to falsification by the upcoming data from the 
LHC and it successors and the large water Cerenkov detectors now being planned to 
further raise the limits on proton decay. This route evades difficulties and confusion 
that had waylaid the fermion fitting program in MSGUTs in the recent past. We have 
shown that the way out is to accept the difficulties as indicative of a particular type 
of sfermion spectra which will modify the d, s quark yukawa couplings significantly 
at the threshold Ms- To our knowledge this is a novel proposal that ties the fermion 
spectrum achievable by a Susy GUT unequivocally to its soft breaking terms. Such 
a program is feasible and plausible precisely because of the tight constraints under 
which the NMSGUT functions and the lack of freedom to add additional fermion 
mass operators as is so abundantly available in SO(10) Susy theories based on low 
dimensional Higgs fields and non-renormalizable operators. 

If the colliders do not belie their promises they may sketch the supersymmetric 
spectra for us. With these spectra in hand the NMSGUT will be have to face both 
edges of the scissor we have shown operative: failure at Mx to fit fermion spectra 
unless the sfermion and gaugino spectrum is such as to push the down fermion Yukawa 
couplings above Ms to values where they match with the otherwise accurate and 
consistent fits we find at Mx- The fits we have found rely on very small 126 couplings 
to boost the Type I seesaw masses while completely killing the Type II masses. 
Interestingly this has the effect of lowering the right handed neutrino masses into 
just the region 10 9 — 10 13 GeV where most Leptogenisis scenarios typically place 
them. Moreover the same fits are consistent with a common GUT scale parameter : 
something that it was feared(from the earliest days of SUSY SO(10) GUTs[5]) would 
definitely suppress Type I masses completely. The no-go mechanism that the work of 



[18] had seemed to have proven against our proposal for ultra small 126 couplings Jab 
has been neatly evaded by the theory in a natural way by accepting the indication of 
the reality of low down type quark masses as a token signal yielded by the model under 
the pressure of the fit. This has the potential to yield information about the sfermion 
masses at Ms and thus about the whole pattern of Susy breaking parameters at 
Mx- Inclusion of non-diagonal corrections and Flavour change constraints dependent 
on all the sfermions will significatly refine the filtration of candidates and push the 
model further towards over determination. When these multiple constraints have 
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been applied to the viable parameter sets emerging from a comprehensive survey of 
the down hill simplex fits that we are now carrying out we expect that the NMSGUM 
will finally truly graduate from a very attractive and persuasive Model to the status 
of a mature Theory (ready to attempt its marriage to gravity ?), if only on the day 
it cannot any longer escape this projected falsifying prosecution and must therefore 
be considered a falsified theory; that day however has yet again been deferred by our 
work presented here. 
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